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Abstract 

A model is proposed to demonstrate that classical general relativity can emerge from loop quantum 
gravity, in a relational description of gravitational field in terms of the coordinates given by matter. 
Local Dirac observables and coherent states are defined to explore physical content of the model. 
Expectation values of commutators between the observables for the coherent states recover the 
four-dimensional diffeomorphism algebra and the large-scale dynamics of the gravitational field 
relative to the matter coordinates. Both results conform with general relativity up to calculable 
corrections near singularities. 

1 Introduction 

Loop quantum gravity [H [6] is a candidate quantum theory of gravity. Its non-perturbative ap- 
proach strictly respect the background independence that is the heart of general relativity. Based on 
Dirac quantization in the Arnowitt-Deser-Misner (ADM) formalism jHEJE], the theory's kinematic 
Hilbert space is rigorously defined and is called knot space. Knot space is spanned by knot states, 
where each knot state is a topological network colored by quantum numbers of gravitational and 
matter fields. The gravitational quantum numbers carried by the edges and nodes of the networks 
give quanta of area and volume in space [H |6]. This quantum geometry successfully describes 
the discretized spatial geometry in Planck scales, and accounts for the Bekenstein-Hawking entropy 
of black holes [10] . 

Because of the background independence, the theory also faces new challenges [HEJ E]: one has 
to solve the intricate Hamiltonian constraint that acts on knot space to obtain a physical Hilbert 
space; and, one has to construct diffeomorphism-invariant local Dirac observables to describe the 
field dynamics. The difficulties make the semi-classical limit of the theory hard to obtain, and it is 
unknown whether general relativity is a semi-clasical limit of loop quantum gravity. 

Following the guidance of loop quantum cosmology j33[ [32], [301 El] , the model in this paper 
obtains its semi-classical limit from knot space by assuming: 1) a modified Hamiltonian constraint 
operator is valid; 2) a group averaging procedure applied to knot space solves the modified Hamilto- 
nian constraint, resulting to the physical Hilbert space HI for the model; 3) the matter back-reactions 
on the gravitational dynamics can be ignored in this context. The matter field operators provide 
spacetime internal coordinates for the gravitational local observables in the model [21j|20j[19|. Ap- 
propriate coherent states in HI are defined to minimize the uncertainty of the gravitational local 
observables, whose expectation values give rise to emergent classical gravitational fields. The sym- 
metry of HI then leads to equations governing the emergent gravitational fields. The equations 
reproduce classical general relativity in the vacuum, up to quantum gravitational corrections and 
matter back reactions. The corrections are expected to be important in small scales or near singular 
regions of the emergent spacetime, and the model also provides means to calculate them in the 
next step. 



1 



2 Loop Quantum Gravity with Matter Fields 



2.1 Knot space 

Loop quantum gravity is based on canonical general relativity in Ashtekar formalism [H O E] . This 
formalism describes gravitational fields in a form similar to that of matter gauge fields. 

Traditionally, canonical general relativity uses spatial metric and extrinsic curvature defined in 
the spatial manifold M as phase space variables [HE1E]. As an alternative, triad fields consisting 
of three orthonormal vector fields {ef (x)}(x = (x, y, z); a = x, y, z; i = 1, 2, 3) can be used in place 
of the spatial metric. In the basis of e"(x) and its inverse e z a (x), the spatial Levi-Civita connection 
and extrinsic curvature take the forms T l a (x) and K l a (x). Ashtekar variables (A % a (x) , Ef (x)) are 
related to (ef (x) , K l a (x)) through a canonical transformation by [21 E]: 

E?(x) = det(e)et(x); ^(x) = F a (x) + jK^x) (2.1) 

where the real number 7 is called the Immirzi parameter. By construction, the fields Ef(x) are 
densitized triad fields, and the fields A l a (x) are 50(3) gauge fields. The variables have the non- 
vanishing Poisson brackets: 

{4(x),£}(y)} = 8^(G/c 3 ) 7 ^(x,y) (2.2) 

where G is Newton's constant. Note that one may replace the 50(3) symmetry group with SU(2) 
in this formalism, since SU(2) and 50(3) share the same Lie algebra. 

Matter fields with a gauge group Q can be included in this formalism |28l l29j [15] . Using the 
triad basis we describe the fermion, scalar and gauge fields by (£f(x), 7ri(x)), (^(x), P\{x)) and 
(AJj(x), Ef(x)). Here i, i and i are respectively (gravitational spin g) SU(2), Q and adjoint Q 
indices. 

In this formalism, the action of general relativity possesses the following local symmetries: 1) 
time re-foliation invariance which introduces the Hamiltonian constraint H(N) = 0; 2) spatial 
diffeomorphism invariance, which introduces the momentum constraint M(V) = 0; 3) local SU (2) 
invariance, which introduces the SU{2) Gauss constraint G(A) = 0; 4) local Q invariance, which 
introduces the Q Gauss constraint G(A) = 0. Each constraint is a functional of the dynamical fields 
and the Lagrangian multipliers N(x), V a (x), A*(x), A : (x) (the bars indicate their non-dynamical 
nature). The first three constraints result from spacetime diffeomorphism symmetry, and consist 
of pure gravitational and matter terms: 

H(N) = H g (N)+H m (N); G(A) = G g (A) + G m (A); M(V) = M g (V) + M m {V) (2.3) 

Denoting n = 87r(G/c 3 ), the three constraints form a closed algebra with structure functionals 
independent of the explicit forms of the matter terms. Denoting [A, A'] to be the SU(2) commutator, 
and setting [V, V'] a = V b d b V' a - V' b d b V a , the algebra is given by 



{G(A),G(A')} = ^G([A,A'}); {G(A), M(V)} = K 7 G(CyA); {G(A), H(N)} = 
{M(V),M(V')} = kjM([V,V'}); {M(V),H(N)} = KjH(CyN) 

[E a d a N, E b d h N' 



{H(N),H(N')} = (M(S) + G(S ■ A)) + ^(1 - 7 2 )G 
where £ v denotes a Lie derivative and 



detE\ 



(2.4) 



rpb Tpai 

S a = (Nd b N' - N'd b N)-f-— (2.5) 

det E 
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In QCD, the generalized electric flux and magnetic holonomy variables are powerful in capturing 
non- perturb at ive degrees of freedom of gluon fields. The Ashtekar formalism enables an analogous 
treatment of gravitational fields. Gravitational holonomy and flux variables are signatures of loop 
quantum gravity [5j 01 |6] , and they capture non-perturbative degrees of freedom of gravitational 
fields. The holonomy variable over an oriented path e C M (the bar indicates that e is embedded in 
the spatial manifold M) gives the parallel transport along the path by the connection fields AUx). 
The flux variable over an oriented surface S C M gives the flux of Ef(x) through the surface. 
Explicitly, we have: 

h®(e)*[A] = [Vex V J_de b Ai(x)r^]f; F<(5) = J_E?dS a (2.6) 

where V denotes path ordering along e, and the 5'C/(2)-valued gravitational holonomy is written 
in the spin (J) matrix representation. The matter variables compatible with the gravitational 
flux and holonomy variables are the following [281 1291 115| . The gauge fields A„(x) are described 
by the Q holonomies h"(e)j in representations (i), and the Ef(x) fields are described by the flux 

variables Fi(5). If v is a generic point in M, the £|(x) fields are described by the irreducible tensors 
9^ d \v) obtained from their Grassmann monomials of degree d, and the spinor momenta vrj(x) are 
described by rj(v) = i0^ d=1 \v)^ . The <^(x) fields are described by h^(v) = exp((fj (v)t^ j) called 
point holonomies in representations (k), and the momenta -P;(x) are described by Pi(v). This new 
set of gravitational and matter variables are collectively called loop variables. The kinematical 
states of loop quantum gravity [5j[4j[6j are called knot states, and they are functionals of the 
configuration Ashtekar variables {^(x), AJj(x), £*(x), ^;(x)} via the corresponding loop variables 

{h®(e),h®@,0W(v),hto(v)}. 

A knot state is given by an SU(2) x Q invariant product of loop variables defined on a graph 
in M, so it solves the Gauss constraints. Further, the knot state is only sensitive to the embedding 
of the graph up to a spatial diffeomorphism \i £ diffM, so it also solves the momentum constraint. 
Define an embedded graph 7 in M to consist of N e smooth oriented paths {ei}, called edges, 
meeting at most at their end points {v n }, called nodes (the bars again indicate that 7 is embedded 
in M). Carrying loop variables, an embedded colored graph T is defined by: 1) an embedded graph 
7; 2) an SU(2) spin representation ji and a Q group representation ] i assigned to each edge; 3) 
generalized SU(2) and Q Clebsch-Gordan coefficients (intertwiners) % n and i n , a point holonomy 
representation k n and a Grassmann monomial degree d n assigned to each node. The assignment of 
jii ii-, in, in, k n and d n gives an SU(2) x Q scalar functional Sf. 

(- N v N v N v N v N e N e \ 

S f [A, A, <f>] = Inv I (g) i n (g) i n (g) 9^ (v n ) (g) (v n ) (g) h,W (e;) (g) (e<) I [A, A, £,</>] 

(. n n n n i i ) 



(2.7) 



where Inv{...} denotes the SU(2) Q invariant contraction. An element fi G diffM drags f to 
T' = and transforms Sp to Spfi = Sfv. Note that for every 7, there is a subgroup diffM, 7 of 
diffM that leaves the graph invariant, maintaining the set of all the edges and their orientations. 
There is also a subgroup TdiffM,^ of diffM, ■y that acts trivially on 7. Thus we have the graph 
symmetry group GM,y = diffM, y/ "I 'diffM, 7- To erase the embedding information in Sf, the group 
averaging operator Fdiff M * s defined as: 



St 



Sr 



N G 



M,7 



A»' 'tdiffM I 'diff M ,7 



M,7 



*[f] 



(2- 
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where Ng m ~ is the number of elements in Gm,^, and [r] is a colored graph obtained from the 
embedded colored graph f by erasing its exact embedding^ The result is a knot state. Such 
a state is determined by a colored graph [T], and is diffM invariant since Srfi • fi = Sr^f] = s [r]- 
The inner products between knot states are given by (generalized) Ashtekar-Lewandowski measure 
HI El [15] , with which the set of all knot states {(srji|} gives an orthonormal basis: 

(s[r']l*[r]) = s [T'},[r] (2-9) 

This basis spans knot space K, the SU(2) x Q and diffu invariant kinematic Hilbert space of loop 
quantum gravity. 

Canonical quantization of loop variables leads to the operators of the form 0(Cl), where Q C M 
may be u, e or S. However, 0(0.) does not preserve K since £1 is not cliff m invariant. To preserve 
cliff m symmetry, we now replace by a dynamical object £1 that assigns O(T) C M to every F, such 
that 0(r) transforms together with f under diffu transformations. Setting {r rep } to contain one 
representative of every [T], one can construct Q by specifying 0(r rep ) on the representatives {r rep }. 
Formally, each dynamical object is a map O : T — > 0,(T) C M satisfying 0,(/jT) = // /j,£l(T) for 
any \x 6 diffu and some fj G Tdiffu,^- Then, we can define a dijffM invariant operator 0(Q) as: 

sm .6(n) = 5 r -6(n(f))p ( i^ (2.10) 

For example, one can construct diffu invariant flux and holonomy operators using a dynamical 
surface S : f — > S(T) C M and a dynamical path e : F — > e(f ) C M: 

S[r] • = 5 r • F(S(f))P di#M ; S[f] • fcCfl( e )f = 5 r • ^)(e(f))fP^ tf (2.11) 

The two operators are respectively differential and multiplicative operators acting on srfj, and their 
SU(2) x Q invariant products give gravitational operators in K. Since Ef(x) determines the spatial 
metric, the spatial area and volume operators are made up of the flux operators [SHUE]. Fi(S(T)) 
as a differential operator acts on the gravitational sector of Sf in a way that satisfies the Leibniz 
rule. Specifically, when S(T) intersects with 7 only at its node v%, we have: 

N v N e _ N v N e 

Fi(S(T)) • (g) i n (g) bf»Xei) = Y, < B > e i .)i{e i .,v x )Ji{e i ,) ■ (g) i n (g) h^^) 
n i ^'kiev n i 

N v N e N v N e 

Ji{e) '®i n ®hM{ei) = ihK^Y, [ S e,e k • h^\e k )T^ k) -8---X ■ T? k) h^(e k )\ (g) (g) ^ (e*) 

n i k n ijkk 

(2.12) 

where c(S, e^) is +1 or —1 when is above or below S in an infinitesimal neighborhood of vi, and 
is if otherwise; i(ei>,vi) is +1 or —1 when v± is the source or target of ey. 



3 The Model 

The remaining constraint to be imposed for a physical Hilbert space is the Hamiltonian constraint. 
Adhering to the polymer-like structure of knot states, the standard Hamiltonian constraint operator 



1 Note that when V = Tp., we have [r] = [r'] 
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H(N) LG:G = Hg® G (N) + H^P G {N) [3 El [15] is quantized from a regularized discrete expression 
approximating H(N). In the discrete expression, the curvature factors in H(N) are approximated 
by holonomies along a certain set of tiny loops. The quantization then leads to H(N) LO - that 
contains holonomy operators based on these loops. Therefore, the action of H{N) LG:G on a knot 
state S[p] involves a change in the graph topology that adds the set of tiny loops to 7. Moreover, since 

a non-constant Lagrangian multiplier N is not diffM invariant, H LG:G (N) does not preserve K in 
general. With a constant N, the action of H LG:G (N) preserves K, but it is intricate and changes the 
topology of the graphs. As a result, constructing a physical Hilbert space annihilated by H LG:G (N) 
is a major challenge, which is currently tackled by both canonical approaches [3j [6l 1221 127j and the 
path integral formalism (6j[7l|8]. In our model, we will modify H LG:G (N) into a graph topology 
preserving operator H(N p ) defined in K. The much simplified setting will allow us to apply group 
averaging method to construct the physical Hilbert space EI of the model, based on certain concrete 
assumptions. 

For each embedded graph 7, we label each of its nodes with an integer n and each of its edges 
connected to the node n by an integer pah|^](n, i). For a given 7, we denote its node n by v%, and the 
oriented path starting from vZ and overlapping exactly with its edge (n,i) by € t ni (fig. la). To each 
pair (e^ j, ■) we assign a minimal oriented closed path i a that lies in 7, containing the outgoing 
path i and incoming path (e^ (fig. la). Using these labels, we define a set of dynamical nodes 
{v m }, satisfying v m (T) = vL m \ such that n(m) is one-to-one. Corresponding to a given {v m }, we 
define a set of dynamical paths {e m j} satisfying e m j(f ) = eL m j )fc ^ and v m (f) = such that 

j(k) is one-to-one. The dynamical closed path e m ,ij is then determined by the outgoing e m .- L and 
incoming e~ ^ dynamical paths. We also define a set of dynamical spatial points p = {pt}, where 
k ranges from to infinity, such that v2. = pk(T) holds for exactly one k value for every T and 
raj^] Notice that there are infinitely many distinct sets of dynamical nodes, paths and spatial points 
satisfying the above. This ambiguity comes from the arbitrariness of identifying the nodes and 
edges between different knot states, due to the absence of a reference background in K. 

The modification from H L< ^ G (N) into H(N p ) - a crucial step for the model - contains two 
elements [41J. First, N is replaced by a diffM invariant lapse function N p (pk), which is a function 
of a set of dynamical spatial points p = {pk}- Recall that the embedded spatial point Pfc(F) 6 M 
transforms under diffM just as f does. That means the classical counterpart of N p (pk) is the field 
iV(x) that transforms as a scalar field. Therefore, in proper semi-classical limits, the operator 
H(N P ) is expected to approximate the diffM invariant H g (N) instead of H g (N). Second, the tiny 
loops (fig. lb) that define the holonomy operators in H LO:G (N) are replaced by the new set of closed 
paths (fig. la) that are contained in the graphs of knot states. By using the corresponding new set 
of holonomy operators, H(N p ) preserves the graph-topology of knot states. Since H(N p ) preserves 
both diffM symmetry and the graph topology, it is an operator in any subspace Kj- C K with 
one specific graph topology T ■ For simplicity, the model's kinematical Hilbert space is set to be 
K%orus C K with the graph topology of a lattice torus Ttorus, which has N v nodes and six edges 
connected to each of the nodes. Restricted to Kj- torua , we have 1 < n < N v and 1 < i < 6 for v n 
and e n> i , and j • will be a square closed path overlapping with exactly four edges in 7 

2 Note that each edge has two labels since it contains two nodes, and the range of i depends on n in 
general 

3 The spatial manifold M contains uncountably many spatial points, but we use only countably infinite 
set {pk} in correspondence to the discrete structure of K. Also, the one-to-one correspondence between v% 
and Pfc(r) is obvious had we taken {pk(T)} for any T to be the set of all points in M. Here we define {pk(T)} 
to be a countably infinite subset of M, while maintaining this natural condition. 
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Figure 1: T/ie Ze/t figure depicts an embedded graph 7, wi/i one 0/ its nodes labeled by n = 1 and t/ie 
/ow edges connected to this node labeled by (1, j). The figures in (a) demonstrate the corresponding 
definitions of Vn, e^- and e^ ni j, while the figures in (6) demonstrate the corresponding definitions 

of el -\ and el . that are used in H LG:G (N). 



The modification applied to the gravitational term of H L( ^ G (N) results to the new gravitational 
term H' g (N p ), which acts on (srfi| G Kj- torus as (set h(e) = h^' 2 '(e)): 



s [f] H' g (N p ) = S f 



p k (f)=v m (t)) ^2 sgn(e m ,<(r),e mj -(r),e m>fc (r)) 

i,j,k=l 



h-\e m ,(t))\ [h{e m>i {t% [^ f) (l), V(f)]] ^ _1 (emj(f))f [h{e m , 3 {t)% [£^(1), Vfo 



where we have 



(3.1) 



^2 N p(Pk\p k (T)=v m (T)) sgn(e mi i(r),e mJ (r),e miA; (r)) 

V m i,j,k=l 

I 



H e m,i,j(T)) - h (e m ,i 7 j(T)))_ [h (e mifc (r)) 



h(e mtk (T))j_,V {t) 



(3.2) 
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Here, the total volume operator Vrf\ is defined as 



%0 = X! 4a ^2 sgn { e ™A T )> e rn,j(r),e m , k (T)) e pqr j p (e mi i(T))j q (e m j(r))j r (e m ,k{T)) 

Vm i,j,k=l 

(3.3) 

The modification applied to the matter term in H L ® G (N) [15] results to the new matter term 
H m (N p ), whose explicit form depends on the matter content. Similar to the pure gravitational 
term, it is also constructed from the loop operators in Kj- torua , and acts on (s^ \ G Kj- torus as: 

s m H' rr XN p ) ^ S T H' m(T) (N p )F dlff 

H'rn(r)( N p) = ^ N p{Pk\ Pk {V)=v n {Y))H'Zyi{^ 

Vn ^ 

h G) (^r))f,h (j) (^,^ 

^(^(r))},p i (^m(f)),p i (^n ) i(f))J 

(3.4) 

where v n ^ denotes the end node of e n ^. 

Finally, the Hamiltonian constraint operator for our model is given by the self-adjoint sum of 
the gravitational and matter terms: 

H(N P ) = Hg(Np) + H m (N p ) = l -{H' g {N p ) + H' m (N p )) + l -{H' g {N p ) + H' m {N p )^ 

Recall that our kinematical Hilbert space K-j- torus is already SU{2) x Q and diffu invariant. To 
obtain the physical Hilbert space of the model, we still need to impose the remaining symmetry 
generated by {exp(iH (N p ))} with arbitrary N p based on arbitrary p. These unitary operators form 
a faithful representation of a group G, that is 



OO 1 

exp(iH(N Pk )) [ = {u(9)} geG (3-5) 



vfc=l J ~K T . 

' torus 



where N Pk is an arbitrary lapse function based on an arbitrary pj~, and ~ K-y torua means that we 
identify two expressions if they give the same operator in Kf torua ■ The physical Hilbert space of the 
model is constructed using group averaging procedure under the assumptions: (1) the existence of 
the left and right invariant measure dg for G; (2) the operator P defined by 

dgU(g) (3.6) 

G 

maps (ip\ G Kj- torua into P|?/>) € K^- . The two conditions hold in minisuperspace models [26" 1 [30 l 
131] . but remain to be proven for the model. Under these assumptions, the inner product between 
any two states = P|^i) and = PIV^} may be defined as 

(*l|*2) = (H^2) = (^l|P|V2) (3.7) 

The physical Hilbert space of the model H C K% oru3 is the space spanned by {P|sjpj)}. By 
construction, HI is invariant under the action of {exp(iH (N p ))} with arbitrary N p , and satisfies 
the modified Hamiltonian constraint. Each element in H is a solution to the quantized Einstein 
equations, and therefore is a quantum state of spacetime with the matter fields. 



7 



3.1 Local Dirac Observables 

The model obtains its local Dirac observables using clocks, spatial coordinates and frames given by 
the matter fields. Given any set of dynamical nodes {v m }, one may build a set of self-adjoint and 
commuting matter operators [H] consisting of scalar operators {(j) (vm), (f> (v m ), 4> 2 {v m ), </> 3 (f m )} 
diagonalized by the knot state basis of Kj- torug , current operators {V}(v m ), Uj(v m )} and conjugate 
current operators {V/ (v m ), Oj{v m )} (I = 1, 2, 3 for the vector currents; 7=1,2 for the spinor cur- 
rents). The operators (f(v m ), 4> 2 {v m ), (fr 3 (v m )) = &(v m ) will serve as spatial coordinate operators, 

{Vf (v m ) , Ufam)} and {V/ (v m ) , (v m )} will serve as spatial frame operators, and <jP{v m ) will be 
the clock operator. 

Classical gravitational fields in the Ashtekar formalism are SU{2) tensors. In the model, the 
SU(2) invariant components of the fields are described relative to the matter spatial frames. Ex- 
plicitly, for any sjjj G Kr torus we define 

S[f] • J(e n ,j)l = S T • V(v n (r))}j(e n:j (T))iF diff 

s [f] ■ h(e n>k ) T j = ■ G(v n>k (f ))ih{e n>k (T)))u(v n (t )) J }p di// (3 ' 8) 

In canonical general relativity with the matter spatial coordinate field c &(x), we can obtain a diffu 
invariant and spatially local variable 0{X) by integrating det(9 ( I ) (x))5( < I>(x) — X)0(x) over M. 
Analogously, the model uses a normalized Gaussian distribution 5 e with finite width e and defines: 

s m 6(X) = 5r^det(A4(« n (f)))5 e ($(S n (f))-X)6(« n (f))P di// 

n 

s^O'iexAx) = det(A$(u n (f ))) det(A*(« n ,<(r )))6 e {$(v n (F )) - X)6 e ($(v n>i (r)) — X — AX) 

n,i 

xd(e n>i (r))P diff 

(3.9) 

for any srji G Kj- torug . where the coordinate volume element operators are given by: 

A4> e- n , l( r)= [*(«n,i(f))-*(« n (r))] 

5rdet(A*(« n (T))) = S f ]T sgn (e n4 (f ), e nj (f), e n , fc (f )) A^.^ • (a^.^ x ^ Sn , k (r)) 

It is crucial that the spatially local operators obtained in (3.9) do not depend on the choices of the 
dummy variables {v m } and {e nj j}. The diffu invariant operators are localized by referring to the 
spatial matter coordinates. 

To specify time using the clock field the model constructs an operator tlx that maps a 
spacetime state |\&) G EI to its spatial slice state TIt\ ^) G Kj- where = T. With to(v) 
denoting an operator acting on v, the operator fir is defined by: 

r_ _ 1 ( 3 - 10 ) 

S[f] n T = S f sym | Yl v^{v n {Y ))8 f -(ct> (v n (T)) - T) \ ¥ diff 

where the symmetrization sym{...} in the ordering of n is applied. The resulting local Dirac 
observables in H are given by: 

6(X, T) = PO(X)il T ; 6'(e XAX ,T) = P6'(e XA x)^T (3.11) 
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In the following, u in (3.10) will serve as an index running over the matter coordinates and frames 
{4>°, <£, f-Vi, f-V 1 , g-Uj, g-U 1 } contracted with the non-zero SU{2) test functions {/*, fy, g % ,gi}- 
The variables v w thus contain information about the momenta of the reference matter fields. 

Finally, applying (3.11) to the fields of concern we obtain the relevant localized observables 
J{ex,AX,T)i, h(e x ,AX,T)j and v u (X,T). 



3.2 Conditions on Matter Coordinates and Frames 

Next, we choose a proper physical state € EI to derive semi-classical limits from the model. 
In order to give sensible descriptions, the localized observables must refer to well-behaved matter 
coordinates and frames. So there are conditions to be imposed on the matter sector of which 
provides the matter coordinates and frames. 

For the clock, we require that any two spatial slices of |^) with cp° = T\ and 4>° = T2 are related 
by causal dynamics, and thus either one of them can be used to reconstruct the spacetime. In the 
model, this condition is imposed as: 

M Tl |tf) + o(h) = ipn T2 |^) + 0{h) = \m) (3.12) 

On each spatial slice IItJ*), we will also impose spatial coordinate conditions. Denoting an 
arbitrary combination of gravitational operators based on {e n ,i} as tp(J(, e n,i)i, h(e n> i)l), we require 
that there exist {w*} and {e* ni \ such that: 

^(j(e;j i ,/ i (e; i jf,/ i t( e ; ! .)pn Tl |^) 

= ^(J(ex n) Ax n j7,^(ex n ,AX„,jJ,/i t (ex„,Ax n , i )5)nT 1 |*) + 0(n) 

with a certain set of values {X n } and {AX m j} satisfying {X m + AJ mj } = {X n }. Here, {f*} 
represents the physical nodes at clock time T\ for an observer using the spatial matter coordinates 
<I>. Naturally, there is also a set of physical spatial points p* = {p* rn \ agreeing with {v^}, such 
that p* m {T) = v^if) for m < N v . Condition (3.13) says that each physical node acquires a 
matter spatial coordinate value, and that the matter frames are physically orthonormal to each 
other at each physical node. The spatial coordinate condition can now be imposed on the map 
( v n-> e ni) ~~ ^ {X n , AX n j). Analogous to the coordinate maps on a torus manifold, the map should 
appear smooth in large scales at most of the physical nodes except at some {v^ b } C {v* t } where 
coordinate singularities occur. We demand |AX nj j| < d for v* t {w^ 6 }, where d is small enough 
that the map appears continuous at large scales. Also, for any m the set {/\X mi i < d} should define 
a parallelepiped in M 3 up to an error of O(d) so the map appears smooth at large scales. Notice 
that once (3.13) is satisfied, the coordinate conditions can be achieved easily through a redefinition 
of coordinates $->$'(!>). 

Using a lapse function N^{ satisfying N^i{p* m ) = J\f{X m ) with an arbitrary function N(X), we 
can apply (3.13) to H g (N#): 

tf s (A^)|*> ^(A^^ 

U g {N) = <p(M(X n ), J(e XnA x n ,)i, h(e XnAXn y 3 ) 

Finally, the well-behaved matter coordinates and frames also lead to the algebraic relations between 
the spatially local operators, when they act on to Lp{J{e* ni )i,h{e* ni )^)ti T ^^ = From now on 
we set (X, AX) G {{X n , AX n ^)}. Also, define a(ey 5 Ay) as the operator localized from d(e nj j) 
which is diagonalized by the basis {(s[f]| £ Kj- torus }. The state (s^\ has eigenvalue +1 if the 
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embedded edge of T overlapping with e n ^{T) has the same orientation as € n ^ (r) and has eigenvalue 
— 1 if otherwise. The algebraic relations are: 

J ] {exAx)iW) = j(exAx)i\<f) + 0(h) 



(3.15) 



J(.e x ,Ax)i,h(e Y ,AY)j W) = 8 X ,Y8AX,AYill r Y(.Ti)jh(eY,AY) I K\ i P) 

-^x,y+Ay^-AX,Ay^p7(T/)^^(ey,Ar)j W) + 0(l 2 h) 
J( e XAx)i,tf(e Y ,AY)j W) =Sx,Y 5 ^XAYilf > l(Ti)jh\e Y ,AY) I K\ l P) 

-8xy+AYS-AXAY il p7( T i) I R h HeY,AY)j\'P) + 0(l 2 p h) 
J (e x , ax) I, J ( e Y, ay) J |</>) = 5x,YbAXAY il ll£iJ K J( e YAY)Ku(eY > AY)\ i P) + 0(l 2 p h) 

These relations enable further calculations after the approximation (3.13) is made. 

3.3 Coherent States and Emergent Fields 

Suppose \^>) satisfies the matter coordinate and frame conditions (3.12) and (3.13), such that the 
local observables { J(exAX, h(e X AX, u u (X, T{)} give meaningful descriptions around the 
moment T\. Since our goal is to obtain semi-classical limits, we impose coherence conditions on 
1^) with respect to the observables: 

J (ex ax, n)^) = (J(ex, AX ,Ti)j>|tf> + 0(ijj); h^XAX^) 1 ^) = (^(e x ,AX,Ti)j)|*) + 0(l 2 p ) 
j> u (X,Ti)\V) = (i> u (X,Ti))\*)+0(h) 

(3.16) 

such that for any two e* j(f ) and ej^ -(f) that share a common node and form a smooth path, the 
expectation values satisfy (recall that |AX nj j| < d for u* { v n b }) : 

{J(ex n AX n ^ T i)i) = (J(ex m AXmj, T i)i) +°{ d ); (h(e Xn AX n ^T 1 ) I J ) = (h(e Xm AX m , J ,T 1 ) I j) + 0(d) 
lfi u {X nt Ty)) = (j>UXn + AX n ^T 1 ))\^)+0(d) 

(3.17) 

Because of the algebraic relations in (3.15), we expect the solutions to (3.16) and (3.17) to exist. 
The explicit construction of flux coherent states in loop quantum gravity is a subtle issue, and it is 
rigorously studied in works such as |35|, I36|. In the following, we will involve only the semi-classical 
properties (3.16) and (3.17) of |\&), which imply that the quantum spacetime has sharply defined 
and approximately continuous values for the local observables at the clock time T\. 

To make contact with classical general relativity, the model maps the expectation values in 
(3.16) to the classical field values through an algorithm, using the matter coordinates as a com- 
mon reference. Suppose {X n } approximatedly occupies a region I 3 C M 3 . Then each physical 
dynamical path e* i is identified with an oriented smooth path e~x„AX ni C I 3 going from X n 
to X n + AJ„ : j. Also, Sx n AX„ t C I 3 is defined to be a rectangular, oriented surface that in- 
tersects ex„AX ni with the same orientation, such that {Bx n ,Ax n ,} defines a cell decomposition 
of I 3 that is dual to the lattice defined by {ex n AX ni }- After the identification, the algorithm 
maps the expectation values {(J(ex,AX, (h(exAX, Ti)j), (u u (X, T{))} to the smooth fields 
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{Ef(X,T 1 ),A{(X,T 1 ),u UJ (X,T 1 )} defined in I 3 . The map obeys the following rulesQ 
/ Ef(Ti)ds a = {J{exAx,Tx)i) + 0(l 2 p ) 

J Sx,AX 

Pexp[ / ^(Ti)(<rj)de 6 ]f = (fe(ex,AX, 2i)f )+0(Zj); ^(X,^) = <z> w (X, Ti)) + 0(ft) 

(3.18) 

The fitting algorithm described above is restricted but non-unique. However, any fit satisfying the 
requirements (3.18) gives a valid correspondence between \^f) and the smooth fields at T\. 



3.4 Semi-Classical Limit of ^ 



Sequentially applying the conditions (3.15), (3.12), (3.16), and (3.18), one finds that (for details 
see HH): 



n-l 



H g {N n ),.... n g (M 3 ), HgWdtHgM) 



{H g (N n ), .... {H g (N 3 ), {Hg^HgiNi)}} ...} 



+ 0(h) + o{d A ) 



(3.19) 



with arbitrary n and Mi- According to (2.4) in the case of an empty matter sector, (3.19) reproduces 
the full (off-shell ) algebra between H g (N), G g (A) and M g (V) in the semi-classical limit of \^>), up 
to the corrections 0(h) + 0(d 4 ). 



■AT 



Moreover, the symmetry H(N P )\*$>) = has two major implications when we choose N p = N$ 
(for details see [IT]): 

(1) With arbitrary n and A/i, the symmetry implies that: 



n-l 



H(N; 



H(Nfi), H(Np),H(N£ 



■No 



IlrJ*) = o 



(3.20) 



One can evaluate the purely gravitational contribution to (3.20) involving only H g (N^{'), which 
can be approximated by % g (Mi) according to (3.19). With this approximation, we may conbine 
(3.19) and (3.20) and find: 

{JE^CiVn), .... {J^CJV-a), {^(JVi), JJ fl (iV-a)}} ...} l^^^r^^i + OCM+OC^+em = (3.21) 

where e m denotes generic matter back reactions given by terms involving H m (N^{ i ). Thus, the 
emergent gravitational fields satisfy the pure gravitational constraints H g (N) = G g (A) = M g (V) = 
up to the corrections 0(h) + 0(d 4 ) + e m . 

(2) In canonical general relativity, the unit lapse function (N(x) = 1) leads to a time foliation 
in which the speed of 4>° is equal to v^o at the spatial slice </>°(x) = T\. Therefore, the time foliation 
using 4>° as the clock is given by iV(x) = l/t^o(x). With M(X) = l/i^o ( X, 7\ ) , the symmetry 
implies: 



H(Nl^) 



0(e x ,AX,T) 







(3.22) 



4 Note that such a map is guaranteed to exist, since we are fitting the smooth fields with infinitely many 
degrees of freedom to the finitely many data points given by the expectation values of the local observables. 
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When 0(ex,AX ,T) is set to be J(ex,AX,T)j and h(ex,AX,T)j, (3.22) gives the clock time dynamics 
of the emergent gravitational fields: 

d 



dT 



Ef(X,T) 

L 

E?(X),[H g (N) + M g (V) + G g (A)] 
+e m + 0(h) + 0(d) 



E?=E? (Ti ) Ai=Al (Ti ) ,N=uJ (Ti ) ,V° =V « (u u (Ti ) ) , A< =A l (i/ w (Ti ) ) 



d_ 



A J b (X,T) 



^(X),[^ g (iV) + M ff (F) + G g (A)] 
+e m + O(ft) + 0(d) 



Ef=E?(T 1 ),Ai=A^(T 1 ),^- 1 (T 1 ),y«=y«K(T 1 )),A»=A»K(T 1 )) 



(3.23) 



Up to the corrections e m + O(h) + 0(d), these are general relativity's equations of motion in the 
Ashtekar formalism, in the gauge {N, V, A} determined by the matter coordinates and frames. 
Specifically, we have: 



V a \ 



iv$=v 



A 1 



f-vj- 



"g-u T 



g-U 1 



I Vg> = V 







f-vj- 



f.yi 



3-Uj 



(3.24) 



g-U 1 



which corresponds to the gauge with the world lines X = const perpendicular to the spatial slice 
at Ti. 

We now look into the correction terms in (3.19), (3.21) and (3.23). For this paper we assume that 
the matter back reactions e m are small and focus on the gravitational corrections, which are: 1) The 
corrections denoted by 0(h) resulting from the uncertainty principle. All the fields are quantum 
mechanical in the model, so there are quantum fluctuations not only in gravitational fields, but 
also in the matter coordinates and frames. Among the corrections, there are terms coming from 
regularization of the inverse triad factor (&etE)~ 1 / 2 appearing in H g by the commutator between 
the total volume and holonomy operators in (3.2). Thus the quantum Hamiltonian constraint, 
which is constructed to be always finite, deviates from the classical Hamiltonian constraint when 
the emergent E fields approaches the scale of l 2 p . Nevertheless, the terms in 0(h) are negligible 
in large E field regions, when the competing classical terms are much greater. 2) The corrections 
denoted by 0(d) resulting from the discretization of space. Remarkably, these terms are of zeroth 
order of h and could still dominate in large scales when the quantum effects are insignificant. Recall 
that d represents the spatial coordinate gap, which is intrinsically finite due to the discretized spatial 
points {p* n }. The effect of the finite value of d comes in two parts. First, the equations governing 
the emergent fields are difference equations that approximate the classical differential equations, 
and therefore it contains corrections of order 0(d). Second, the use of holonomies instead of the A 
fields as configuration variables introduces corrections of higher orders of \A\d. In the regularized 
classical Hamiltonian constraint, the nonlinear terms in the holonomies contribute to errors, which 
vanish in the limit of the paths shrinking to zero length. In the quantum theory, space is discrete 
and d is finite, and thus the corrections are intrinsically finite. Suppressed by small d value when 
A fields are bounded, the holonomy corrections may become important when A fields become 
singular. Overall, we see that all the corrections denoted by 0(d) are negligible when the emergent 
gravitational fields are nonsingular and vary nicely across a coordinate gap d. 
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When all these corrections are suppressed, general relativity emerges from the semi classical 
limit of the state £l. It should be emphasized that the model shares the kinematics of loop 
quantum gravity, and the modified Hamiltonian constraint H{N p ) preserves most of the original 
features. As a result, the key features of loop quantum gravity- quantum geometry, inverse triad 
corrections and holonomy corrections- are all faithfully carried by the model. Moreover, the model 
provides a set-up to explicitly calculate the correction terms for the dynamics of emergent gravi- 
tational fields, including matter back reactions. It is of great interest to see how these corrections 
behave near the initial singularity, or near a black hole. In mini or midisuperspace quantum cos- 
mology incorporating the key features of the loop quantum gravity- loop quantum cosmology- the 
quantum and holonomy corrections have been extensively studied. Among the remarkable results 
of these models is that the initial singularity can be replaced with a well-behaved bouncing of the 
scale factor, preceded by a contracting universe |_33j|32j and followed by a built-in slow-roll infla- 
tionary phase [37J [38 1 [39J . These predictions may be testable, and give distinguishable signals in 
the spectrum of cosmic microwave background radiation [37 [38j[39j. It is then important to ask 
whether loop quantum cosmology emerges from a certain symmetrical semi-classical limit of loop 
quantum gravity, and what additional details the full theory would provide regarding the predicted 
signals. Explicitly evaluating the correction terms in our model could provide answers to these 
questions. A concrete, affirmative result for the scale- factor bouncing has been obtained in [42J. 
Along with many other calculable physical predictions, the model serves as a promising testing 
ground for ideas from loop quantum gravity. 

4 Conclusion 

The model proposed here has the same kinematics of loop quantum gravity coupled to matter fields, 
which is given by knot states. To bypass the difficulties faced by the standard approach, the model 
uses the modified graph-preserving Hamiltonian constraint operator. This modified constraint 
operator enables the construction of the model's physical Hilbert space through a group averaging 
procedure, based on concrete assumptions. Strictly respecting background independence, the model 
utilizes its matter sector to provide the coordinates and frames to describe the local gravitational 
observables. The coherent state |\&) for the local gravitational observables give rise to classical 
gravitational fields which obey vacuum general relativity up to the matter back reactions and 
the quantum gravitational corrections. The corrections have clear interpretations, and the model 
provides a set-up to explicitly calculate them. It would be a valuable next step to see how these 
corrections behave near the initial singularity, and near the singularity of a black hole. 
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